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Abstract 

We study the phase diagram of the three-state Potts model on a 
triangular lattice with general interactions (ferro/antiferromagnetic) 
between nearest neighbor spins. When the interactions along two 
lattice-vector directions are antiferromagnetic and infinitely strong, 
this model becomes equivalent to a six-vertex model and exhibits 
a first-order (KDP) transition from an ordered phase into a criti- 
cal phase. Comparing the excitations occurred by relaxing the re- 
striction of infinite-strength interactions and those in the eight-vertex 
model, we analytically obtain the critical index for those excitations 
and demonstrate the existence of a critical phase for the case of finite 
antiferromagnetic interactions in two directions and ferromagnetic in- 
teractions in the other direction. When the interactions are antiferro- 
magnetic in all three directions, Monte Carlo simulations show that a 
first-order line emerges from the KDP point and separates completely 
an ordered phase and a disordered phase. Along the special line where 
all three antiferromagnetic interactions have the same strength, the 
cell-spin analysis reveals that the symmetry of the ground states is 
dual to the symmetry of the n = 3 ferromagnetic cubic model which 
is known to exhibit a first-order phase transition. 

1 Introduction 



The ferromagnetic three-state Potts model has been studied extensively. In 
two dimensions, its critical properties which are independent of underlying 
lattices are known exactly by the extended scaling^ || and/or the confor- 
mal invariance . When interactions between neighboring spins become an- 
tiferromagnetic, the critical properties vary with the structure of underlying 
lattices. The symmetry of the antiferromagnetic ground states is constrained 
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by the structure of the underlying lattice. For example, the antiferromag- 
netic three-state Potts model on a square lattice is disordered at all tem- 
peratures, but on a triangular lattice a first-order phase transition appears 
at a finite temperature. By adding ferromagnetic next-nearest-neighbor in- 
teractions, the antiferromagnetic three-state Potts model on a square lat- 
tice exhibits a sequence of two Kosterlitz-Thouless (KT) transitions Qj. This 
model possesses a ground-state symmetry similar to the ferromagnetic six- 
state clock model||. In the case of mixed-type interactions, i.e. ferromagnetic 
in one direction and antiferromagnetic in the other direction, Monte Carlo 
simulations and transfer matrix calculations 0] on a square lattice indicate 
that there is a KT-like infinite-order phase transition from a massless low- 
temperature phase into a disordered phase with an essential singularity at 
a finite temperature. In this paper we examine the phase diagram of the 
three-state Potts model with antiferromagnetic and mixed-type interactions 
on a triangular lattice. 

The three-state Potts model on a square lattice can be mapped to a 27- 
vertex model on its dual lattice Q (section 2). The triangular-lattice model 
with nearest-neighbor interactions can be viewed as the square-lattice model 
with nearest- and next-nearest-neighbor interactions by distorting the lattice 
properly. We map this model onto a 27-vertex model on its dual lattice. 
When the interactions along two lattice-vector directions are antiferromag- 
netic and infinitely strong, only six vertex configurations survive. From the 
exact solution of the six- vertex model ||, we find a first-order (KDP) tran- 
sition from an ordered phase into a critical phase. By the stability analysis 
similar to the work by den Nijs et aZ@] for the square-lattice model, we 
demonstrate the existence of a critical phase for the case of finite antiferro- 
magnetic interactions in two directions and ferromagnetic interactions in the 
other direction. It implies that there is a KT-type transition from a critical 
phase into a disordered phase in the case of mixed-type interactions. 

The three-state Potts model with the isotropic antiferromagnetic interac- 
tions has been studied previously by the real-space renormalizationH], series 



expansions [ID], and Monte Carlo simulations |]1 If. It has been shown to ex- 
hibit a strong first-order transition. In section 3, its ground-state symmetry 
is investigated. We construct the cell-spin hamiltonian by calculating domain 
wall energies and show that this model can be renormalized to the dual model 



of the n = 3 cubic model [0. The calculated values of coupling constants of 



this cubic model guarantee that the phase transition of our isotropic model 
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is of first order. 

When the interactions are anisotropic, the chirality in domain wall ener- 
gies appear. In section 4, we perform Monte Carlo simulations to understand 
the role of the chirality in this model. We find that a first-order line emerges 
from the KDP point and separates completely the antiferromagnetic ordered 
phase and the disordered phase. The chirality neither changes the nature 
of the phase transition, nor gives rise to any other phase transition in the 
antiferromagnetic region. Pleliminary results for the antiferromagnetic case 



have been published separately elsewhere [J3 



We conclude in section 5 with a brief summary. 



2 Mapping to the vertex model and stability 
analysis 

Consider the three-state Potts model on a square lattice with nearest-neighbor 
and next-nearest-neighbor interactions. The hamiltonian of the model, H, is 
given as 

-H=Y, KS om + £ LicW, + E £a«W (!) 

where < .. > and (..) denote nearest neighbors and next-nearest neighbors 
respectively (see Fig.l). Oi is the Potts spin at site % which takes the values of 
0,1,2 and 5 is the Kronecker delta function. The three-state Potts model on 
a triangular lattice with nearest-neighbor interactions only can be obtained 
by taking either L\ or L 2 to be zero. Here we take L 2 = 0. 

A three-to-one mapping to a 27-vertex model on the dual lattice is ob- 
tained by assigning arrows or zeros to the bonds of the dual lattice. As one 
goes around a dual lattice site clockwise, an outgoing (incoming) arrow is 
assigned on the encountered bond if the value of the Potts spin is increased 
(decreased) by one with modulo 3 going across the bond. If the value of 
the Potts spin is unchanged, we assign a zero on the encountered bond (see 
Fig.l). Bolzmann weights in terms of u = exp(K) and v = exp(Li) and the 
number of vertices are shown in Fig. 2. In the limit K — ► — oo (u — 0), only 
six vertices are left with nonvanishing Boltzmann weights. By normalizing 
Boltzmann weights with respect to the last pair of unpolarized vertices, we 
find a = 1/v and 6 = 1 where a and b are the Boltzmann weights of the first 
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two pairs of vertices. The single parameter A of the six- vertex model [EL [14 
is given as 

A - - (- - - —\ — (2) 



From the exact solution of the six- vertex modelH [14J, we find that there is 
a first-order KDP transition at A = 1 (v — 1/2) from an antiferromagnetic 
ordered phase (v < 1/2) into a critical phase (v > 1/2) (see Fig. 3). 

With finite K (u > 0), the vertices with zeros on the bonds can appear. 
Following the analysis of den Nijs et a/|3J for the antiferromagnetic three- 
state Potts model on a square lattice, the most important vortex excitations 
which drive the critical phase into the disordered phase are the bound pairs 
of vortex states (Figs. 2(e) and 2(f)) which have vorticity of ±6 (Fig.4). 
The scaling dimension of these excitations can be obtained by using the 
well-known relation between scaling dimensions for excitations of different 



vorticities such as 15 



m x 2 



•^m [ ) ■^ni (3) 

where x m is the scaling dimension for excitations of vorticity m and zero spin- 
wave excitation index. The scaling dimension is known exactly from the 



Baxter's solution of the eight- vertex model|T6|. Thus the scaling dimension 
for excitations of vorticity 6 is given as 

x e = (V) %8V, (4) 



where x&v = 2 — y$v = 2 — |cos _1 (— A). The critical exponent y$ = 2 — 
Xq becomes negative when v > v* = 1/(2 sin(y|)). So these excitations 
are irrelevant in this region and the critical phase persists for a small but 
finite value of u. For v < v*, they are relevant with respect to the u = 
line. In this region, these two bound pairs can be dissociated and the 
system becomes disordered for any finite value of u. Thus we can draw 
the phase diagram for small u in the axis of v (see Fig. 5). As v increases 
from the antiferromagnetically ordered phase, a first-order transition into 
the disordered phase is expected near v ~ 0.5 and subsequently a continuous 
KT-type transition into the critical phase near v ~ 2.879. Notice that the 
bound pairs of vortex states in Fig.4 are always confined for small u because 
a string of zeros are generated by pulling the bound pairs apart. 
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3 Ground state symmetry and cell-spin hamil- 
tonians 



First we study the ground state symmetry of the three-state Potts model 
with antiferromagnetic nearest-neighbor interactions and ferromagnetic next- 
nearest-neighbor interactions on a square lattice^, Its hamiltonian is 
given in eq.(|l|) with K < and L = L\ = L 2 > 0. This model has six 
equivalent ground states. The unit cells of these ground states are shown 
in Fig. 6. At zero temperature, the system becomes a periodic array of one 
of these unit cells. As the temperature goes higher, the system becomes a 
mixture of these unit cells and the domain walls between different unit cells 
appear. The excitation energies per unit length for these domain walls, when 
they are straight and do not meander, are easy to determine 

E hi+1 = L, E iti+2 = 2L-K/2, and E iti+2 = 2L-K, (5) 

where i is an integer of modulo 6 and Ey is the excitation energy of the 
domain wall between two ground states, % and j. There is no chirality in 
domain wall energies; Ey = Eji. We observe from eq.([5D that the domain 
wall energy Eij depends on \i — j\ only and is a periodic function of i and 
j with periodicity of 6. This is a symmetry of the six-state clock model 
where the domain wall energies depend only on the angle between states 
(Fig. 7). If we assign a cell spin a (a — 1, • • • , 6) for each unit cell, the cell- 
spin hamiltonian reduces to the ordinary six-state ferromagnetic clock model 
hamiltonian 

-H = - J2 JIl-coA* (6) 

<i,j> 

with L = —K/2 = J/2. The six-state clock model is known to exhibit a 



sequence of two KT transitions [IT]]- This explains why there exists a critical 
fan in the antiferromagnetic three-state Potts model on a square lattice. 

Now we consider the three-state Potts model with antiferromagnetic 
nearest-neighbor interactions on a triangular lattice. The coupling constants 
along three different lattice- vector directions are denoted by Ki (i = 1, 2, 3). 
All Ki§ are negative. Similar to the above square-lattice model, there are 
six equivalent ground states; three up-states Ui and three down-states A: 
(i = 1,2,3). Their unit cells are shown in Fig. 8. We say that the up-states 
have a positive helicity and the down-states a negative helicity. Domain wall 
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energies between these ground states depend on their directions and also on 
the chirality. Consider the domain wall between two up-states U\ and U2 
in the direction 1 (Fig.9). When U\ is left and U 2 is right to the domain 
wall, the excitation energy per unit length of this domain wall is —K 3 . If the 
ground states are interchanged, the domain wall energy becomes —K 2 . So 
when K 2 7^ K 3 , there is a chirality in the domain wall energies. In general, 
one can find that the domain wall energy between two up-states C/j and Uj 
in the direction 1 is E\ — —K-ih-xj — When % = j — 1 (J + 1), we 

call that the domain wall has a positive (negative) chirality. Similarly, the 
domain wall energies between two up-states in the other directions can be 
obtained easily and the result is 

Ef{UU) = (7) 

where the subscript i denotes the direction of the domain wall, the superscript 
± the chirality, and UU in the parenthesis represents the domain wall energy 
between up-states. Repeating the same analysis on the domain walls between 
down-states and also between up- and down-states, we find 

Ef{DD) = -K i±1 , 

Ef{UD) = -l^+. + K^). (8) 

Notice that the domain wall energies between up- and down-states do not 
depend on the chirality. 

The symmetry structure of the six ground states is drawn in Fig. 10. There 
is a ferromagnetic chiral (or helical) three-state Potts model symmetry [|18|1 
in each triangle. And these triangles are linked by the symmetry of a ferro- 
magnetic nonchiral Ising model. When K\ = K 2 = K 3 (isotropic case), the 
chirality disappears in both triangles. Even though the number of ground 
states is the same as in the square-lattice model discussed previously, the 
symmetry between the ground states is completely different from each other. 

First consider the K\ = K 2 = — 00 (u = 0) limit. This is the six- vertex 
model limit (see section 2). In this limit, only four walls can survive and 
their energies are 

Et{UU) = E 2 {UU) = E^(DD) = E 2 (DD) = -K 3 . (9) 

Domain walls in the direction 3 are not allowed and the up-states and down- 
states cannot coexist. There will be no isolated loop excitations of domain 
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walls because there exist no pairs of domain walls in the same direction with 
different chirality. Thus any excitation of this model can be represented 
by the zig-zag lines of the domain walls of types E^(UU) and E^iUU), 
or E^(DD) and E^iDD) (Fig. 11). The domain wall energy is given by 
v = exp(K 3 ) and there is no energy cost at the crossing of the walls. By 
mapping to the vertex model with Bolztmann factors normalized with respect 
to the last pairs of the vertices, one can recover the six-vertex model with 
a = 1/v and b = 1 as expected (Fig. 12). So in the u = limit, the cell-spin 
approach produces the exact result. 

For the isotropic model (u = v = exp(K)), there is no chirality. Assign 
two types of cell spins, t and s (t = 1,2,3 and s = 1,2), for each unit cell. 
The s = 1 state represents the up-states and the s = 2 state the down-states. 
Each of three states inside the up- or down-states is represented by the spin 
t. Then the cell-spin hamiltonian can be written as 

- H = E \-K8 SiSj 8 titj + ^5 SiS . + \k] , (10) 

<i,3> 1 6 6} 



which is exactly the same as the hamiltonian of the so-called (q s , q t ) model [19 
with q s = 2 and q t = 3. For q t — 2 it is known as the cubic model |T2|]. Nature 
of phase transitions does not depend on the underlying lattice structure for 
ferromagnetic models like the above cell-spin model (-K > 0). So we study 
the above model on a square lattice which has been investigated in details. 
The duality relation between the (q s , qt) model and the (qt, q s ) model is known 
on a square lattice|19[. After dropping the constant term in eq. (|T0D , one can 



find the dual hamiltonian 

-H D = [D6 titj 6 SiSj + J5 titj ] , (11) 

<i,3> 

where 

6 



exp(D) = 1 



exp(-2K/3) + 2exp(K/3) - 3 : 



This is the n = 3 cubic model hamiltonian which is known to exhibit a first- 
order phase transition for J + D/2 > and a continuous phase transition 
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otherwise [T^] . For our model, we can prove from eq.(|l2|) that D > and 
J + D/2 > for any value of K. Therefore the cell-spin analysis for the 
antiferromagnetic isotropic three-state Potts model on a triangular lattice 
shows that there must be a first-order transition rather than a continuous 
transition and the symmetry of the ground states is dual to the symmetry of 
the n = 3 cubic model. Nature of the transition is consistent with the Monte 
Carlo results [ |TT|1 . 

On the other side of the phase diagram (K < 0,L > 0), there exists a 
critical phase. There are infinitely many ground states in the thermodynamic 
limit. In the ground states, Potts spins are ordered completely along the di- 
rection 3 and nearest-neighbor spins in the other directions should be different 
(Fig. 13). This model may exhibit the same critical behavior as the square- 
lattice model with mixed-type nearest-neighbor interactions only; antiferro- 
magnetic in one direction and ferromagnetic in the other direction || [7|, ^0]. 
The number of ground states grows exponentially with the linear system size 
N; ric = 2 N + 2(— 1)^ which is much smaller than the square-lattice model 
with antiferromagnetic interactions in both directions. So this model may 
not be disordered at all temperatures in contrast to the square-lattice aniti- 
ferromagnetic model. Our stability analysis in the previous section suggests 
that the transition into the critical phase is of KT type, which is consistent 
with numerical results by Monte Carlo simulations and transfer matrix 
calculations for the square-lattice model with mixed- type interactions. 



4 Monte Carlo simulations and phase dia- 
gram 

Consider the antiferromagnetic region of the phase diagram (0 < u, v < 1). 
Along the u = v line (isotropic case), it is shown in the previous section that 
the antiferromagnetic three-state Potts model on a triangular lattice can be 
renormalized, by the cell-spin approximation, to the dual model of the n = 3 
cubic model and then should exhibit a first-order phase transition. 

For the anisotropic model, the chirality between ground states appear. 
There has been some interests in the role of the chirality in the ordinary fer- 



romagnetic Potts models^, 21|. Even though the chiral operator is relevant 



for the three-state Potts model, it does not introduce a new independent ex- 
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ponent, i.e. xch = %t + 1 where xqh and xt are the chiral and temperature 
scaling dimension respectively^]. So it is believed that the scaling behav- 



ior does not change in the presence of the chirality for the three-state Potts 



model. This has been shown analytically for the hard hexagon model [ 23 1 and 



numerically for the chiral three-state Potts model ]18) and the triangular Ising 



lattice gas[|4[]. The ground-state structure of our anisotropic model is much 
more complicated than that of the chiral three-state models. As explained in 
the previous section, it has two chiral three-state model symmetries linked 
by the Ising symmetry. So it may be quite interesting to find out what kind 
of role the chirality plays in this model. We use Monte Carlo simulations to 
investigate the phase diagram of this model. 

We run conventional heat-bath Monte Carlo simulations on a 60 x 60 
triangular lattice along the u = v, u 2 = v , u = v 2 , and u 3 = v lines. 
Typically a few 10 4 Monte Carlo steps per spin (MCS) are performed for a 
given value of u and v. We measure the antiferromagnetic order parameter 



m^F |nj and the energy density e. The order parameter tuaf is defined as 

where is the number of spins of state a in the sublattice X and N t 
is the total number of spins. Unfortunately, along all four lines, we find a 
very strong first-order phase transition from an antiferromagnetically ordered 
phase into a disordered phase. In Table I, we list the values of the coupling 
constant u at the first-order transitions and the jump of the order parameter 
and the energy density. The five first-order transition points (including the 
KDP point) can be connected by a smooth line which starts from the KDP 
point (u,v) = (0, |) and apparently ends at the point (0.28,0) (Fig.14). 
As the anisotropy increases, the magnitude of the order-parameter jump 
becomes bigger. Our results imply that the strong first-order transition due 
to the cubic nature of this model preempts all possible continuous transitions 
in the whole antiferromagnetic region. 

We also run Monte Carlo simulations in the region of mixed-type interac- 
tions where a critical phase is expected. The specific heat is measured along 
the u = 1/v line. We find a characteristic of the KT transition in the shape of 
the specific heat, which has a very broad bump well outside of the expected 
critical phase. This is consistent with our analytical result which suggests 
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the existence of KT transitions in section 2 and also with some numerical 
results for the square-lattice model with mixed-type interactions || 0, p0[| . 



5 Summary 

We investigate the phase diagram of the three-state Potts model on a trian- 
gular lattice with ferro/antiferromagnetic nearest neighbor interactions. In 
the limit of infinite antiferromagnetic interactions along two lattice-vector 
directions, this model maps onto the six-vertex model and we find the KDP 
first-order transition from an antiferromagnetic ordered phase into a criti- 
cal phase. By the stability analysis, we demonstrate the existence of the 
KT transition from a critical phase into a disordered phase for the case 
of finite antiferromagnetic interactions in two directions and ferromagnetic 
interactions in the other direction. The critical parameter in the limit of 
infinite antiferromagnetic interactions is analytically obtained. In this case 
the ground states are equivalent to those of the three-state Potts model on 
a square lattice with mixed-type interactions. Our result implies that this 
square-lattice model should have the same kind of the KT transition, which 



confirms previous numerical results [20 



When the interactions are antiferromagnetic in all three directions, we 
find six ground states. These ground states have two chiral three-state model 
symmetries linked by the Ising symmetry. For the isotropic model, the chi- 
rality disappears and the cell-spin analysis reveals that the isotropic model 
can be renormalized to the dual model of the n = 3 ferromagnetic cubic 
model with coupling constants guaranteeing the first-order phase transition. 
Monte Carlo simulations for the anisotropic model shows that a first-order 
line emerges from the KDP point and separates completely the antiferro- 
magnetic ordered phase and the disordered phase, which indicates that the 
chirality is not relevant in this model. 
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Table Caption 



Table 1 : Numerical values of the coupling constant Ut, the order-parameter 
jump Aitlaf, and the energy-density jump Ae at the first-order tran- 
sitions. Numbers in parentheses represent the errors in the last digits. 





u = V 


u 2 = V 


u = v 2 


u 3 = V 




0.205(1) 


0.263(2) 


0.120(2) 


0.278(3) 


Attiaf 


0.70(3) 


0.81(3) 


0.78(2) 


0.88(2) 


Ae 


0.17(2) 


0.20(2) 


0.21(1) 


0.24(1) 
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Figure Captions 



Fig.l : A typical configuration of Potts spin on a square lattice (thick line) 
and its corresponding vertex configuration on its dual lattice (thin line) . 
Coupling constants of nearest-neighbor (K) and next-nearest-neighbor 
interactions (L 1; L 2 ) between Potts spins are shown. 

Fig. 2 : The 27-vertex respresentation on a square lattice for the three-state 
Potts model on a triangular lattice: vertex types, their Boltzmann 
factors, and the number of each vertex type. 

Fig. 3 : Phase diagram of the six- vertex model. The dashed line corresponds 
to our model at u — 0. 

Fig. 4 : Bound pairs of vortex states in Figs. 2(e) and 2(f). 

Fig. 5 : Phase diagram of our model for small u. 

Fig. 6 : Unit cells of the six ground states of the square-lattice model. 

Fig. 7 : Ground-state symmetry of the ferromagnetic six-state clock model. 

Fig. 8 : Unit cells of the six ground states of the triangular-lattice model: 
Ui (up-state) has a positive helicity and Di (down-state) a negative 
helicity. 

Fig. 9 : The domain wall between Ui and Uj ground states in the direction 
1. Here i' = i + 1, i" = i - 1, f = j + 1, j" = j - 1. 

Fig. 10 : Ground-state symmetry of the triangular- lattice model. 

Fig. 11 : Domain walls (dashed lines) on a distorted triangular lattice at 
u — 0. The thick bonds are unhappy (ferromagnetic) bonds which 
domain walls go across. Notice that there are no unhappy bonds where 
two domain walls cross. 

Fig. 12 : Six vertex types of domain walls (dashed lines) with unnormalized 
Boltzmann factors. Solid lines represent no domain walls. 
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Fig. 13 : One of the infinitely many ground states for the model with mixed- 
type interactions. 

Fig. 14 : Phase diagram of the antiferromagnetic three-state Potts model on 
a triangular lattice. 
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